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Abstract 

We present the analytic solution of the self-organized critical (SOC) forest-fire 

model in one dimension proving SOC in systems without conservation laws 

by analytic means. Under the condition that the system is in the steady state 

and very close to the critical point, we calculate the probability that a string of 

n neighboring sites is occupied by a given configuration of trees. The critical 

exponent describing the size distribution of forest clusters is exactly r = 2 

and does not change under certain changes of the model rules. Computer 

simulations confirm the analytic results. 
PACS numbers: 05.40.+j, 05.70.Jk, 05.45.+b 
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The concept of self-organized criticality (SOC) [p]] has attracted much attention during 
the last few years since it might explain the origin of fractal structures and of 1/ /-noise in 
nature. While the sandpile model [[J, which is the prototype for SOC, is well understood 
and has also been investigated analytically 0, the investigation of other SOC systems is 
mainly restricted to computer simulations In particular the mechanism leading to 

SOC in non conservative systems is barely understood. So far, there exists no proof and 
even has been questioned that such non conservative systems can become critical. 

A recent paper introduced a critical forest-fire model which is non conservative and 
showed that a double separation of timescales leads to SOC 0. A series of computer 
simulations confirm the criticality of the model |8|-|lT|, and a scaling theory yields relations 



between various critical exponents [flQ|Jll[] . An analytic proof for the criticality of this model, 
however, has not been given so far. || contains the mean-field theory of the forest-fire model, 
which seems to agree with simulations in high dimensions. Critical exponents for the one- 



dimensional forest-fire model have been derived by another mean-field approximation in |T2 
but do not agree with simulations. 

In this paper, we present the analytic solution of the one- dimensional forest-fire model, 
thus proving for the first time that a non conservative system can indeed show SOC. The 
analytic calculation shows also that the value of the critical exponent which characterizes 
the size distribution of forest clusters does not change under certain changes of the model 
rules. The values of the critical exponents have already been given in [[/J by an argument 
which was rather intuitive and failed in higher dimensions. 

The forest-fire model |7j is defined on a (^-dimensional hypercubic lattice with L d sites. 
Each site is occupied either by a tree or a burning tree, or it is empty. The state of the 
system is parallely updated by the following rules: 

(i) A burning tree becomes an empty site. 

(ii) A green tree becomes a burning tree if at least one of its nearest neighbors is burning. 

(iii) A green tree becomes a burning tree with probability / <C 1 if no neighbor is burning. 



(iv) At an empty site a tree grows with probability p. 



A forest-fire model without rule (iii) was introduced earlier |13| 



After a transition period, the model assumes a steady state with a constant mean forest 
density p. Then the mean number of growing trees equals the mean number of burning 
trees. During one timestep, there are on an average fpL d lightning strokes in the system, 
and p(l — p)L d trees grow. Consequently, the mean number of trees destroyed by a lightning 
stroke is 

s={flp)-\l-p)lp. (1) 

When the tree growth parameter p is small enough for a given value f/p, forest clusters 
that are struck by lightning burn down before new trees grow at their edge. In this case the 
dynamics of the system are completely determined by the ratio f/p. In the limit f/p—*Q, 
s diverges, and the system approaches a critical point, where the size distribution of forest 
clusters n(s) obeys a power law. Since any finite system is never exactly at the critical 
point, the power law breaks down at a cluster size s max oc (/ /p)~ x , and is given by n(s) oc 
s~ T C(s/s max ). The cutoff function C is C(x) — 1 for x < 1 and C(x) — > for large x. 
The critical behavior in the forest-fire model can be observed independently of the initial 
conditions and over a wide range of parameter values, as long as / <C p with sufficiently 
small p, i.e. as long as time scales are separated. Then the forest-fire model is self- organized 
critical |J. 

We now turn to the one-dimensional case of this model and investigate the properties 
of its steady state. To avoid finite-size effects, we always assume that the system size L 
is L ^> p/f. The condition of time-scale separation reads in one dimension p <C f/p. 
Consequently there are never two lightning strokes or two growing trees at the same time 
within a distance ^ p/f. Let n(s) be the mean number of forest clusters of s trees, divided 
by the number of sites L. In M, we already derived the following relations 

oo 

Y^n(s) = (l-p-pf/p)/2, 



oo 



^sn(s) = p, 



(2) 



8 = 1 



OO 



Y / s 2 n(s) = (l-p) P /f- 



The first relation follows from the condition that the total number of forest clusters is 
constant in the steady state. 

Before continuing, we need some definitions. Each site can be in two possible states 
which we denote (empty site) and 1 (tree). Let i a £ {0, 1} be the state of site a and let 
P n (ii, . . . ,i n ) be the probability that an arbitrarily chosen string of n sites is in the state 
(ii, . . . ,i n ). We also define the conditional probabilities f n (i 1 , . . . , i n ) that the configuration 
(ii, . . . , i n ) is ignited by fire during one timestep. f n (ii, . . . , i n ) is the sum of the probability 
that a tree belonging to the configuration (ii, . . . ,i n ) is struck by lightning, and, if %\ — 1 
or i n — 1, the probability that fire enters the configuration from outside. In the following, 
we calculate the probabilities P n and f n neglecting terms which are smaller by a factor f/p 
than the leading terms. The results become more and more exact when the critical point is 
approached. 

The probabilities P\{i) can be expressed in terms of the mean forest density: Pi(l) = p 
and .Pi(O) = 1 — p. At any given site, growth and burning of trees occur equally often, and 
therefore 



/i(l) is the sum of the probability that a tree is ignited by a neighbor and the probability / 
that it is struck by lightning. The latter is smaller by a factor f/p than /i(l) and therefore 
negligeable. 

In addition to symmetry relations, the probabilities ^2(^1,^2) satisfy P2(00) + i-2(10) = 
1 - p and P 2 (10) + P 2 (H) = p. The probabilities / 2 (M 2 ) are related by / 2 (10)P 2 (10) + 
/ 2 (11)P 2 (11) = /i(l)Pi(l) = p(l — p) and, since the system is in the steady state, by 
/ 2 (11)P 2 (11) = 2pP 2 (10) and (p + / 2 (10))P 2 (10) = pP 2 (00). When a tree is ignited, the fire 
comes either from the right neighbor or from the left neighbor. This leads to p(l — p) = 



h(l)=p(l-p)/p. 



(3) 
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/ 2 (11)P 2 (11) and / 2 (10) = 0, and finally to 

P 2 (00) = P 2 (10) = (1 - p)/2, P 2 (ll) = (3p - l)/2. (4) 

Remember that terms of order f/p have been neglected. In the limit f/p — > 0, i.e. at the 
critical point, Eq. (|]) becomes exact. Before we consider the case n > 3, we would like to 
comment the results Eq. (|j). 

(i) P 2 (10) = J2^Li n ( s ) ( see (2)) since the configuration 10 represents the right-hand edge 
of a forest cluster. 

(ii) The exact result for / 2 (10) is / 2 (10)P 2 (10) = fp since each time lightning strikes the 
system a right-hand edge of a forest cluster burns down. / 2 (10) is smaller by a factor 
f/p than / 2 (11) and therefore negligeable since only a portion oc f/p of all clusters 
burn down before trees grow at their edge. Below we will see that the approximation 
/ 2 (10) = leads to f n {l . . . 10) = for all n > 2 simplifying considerably the calcula- 
tions. f n (l . . . 10) = cannot be true for large values of n since forest clusters of size 
p/f are struck by lightning with a large probability, which leads to /^(l ... 10) ^ 0. 
Our subsequent results will therefore only be valid for n <ti p/f. 

We are now 

able to calculate the probabilities P n . From f n {l . . . l)P n (l . . . 1) + f n {l . . . 10)P n (l . . . 10) = 
/ n _a(l . . . l)P n _x(l ... 1) and /„(1 . . . 10)P n (l ... 10) < f n ^(l . . . 10^(1 ... 10) we obtain 
by recursion f n (l . . . l)P n (l . . . 1) = p(l — p) and f n (l . . . 10) = 0. This means that a string 
of n sites is covered by a completely dense forest when it catches fire. This dense forest 
belongs to a large forest cluster which has a mean size s > n at the moment when it is 
struck by lightning. As long as not all sites are occupied by trees, the dynamics of our 
string are completely determined by tree growth. All configurations which contain the same 
number of trees therefore have the same probability. Let P n {m) be the probability that the 
string is occupied by m trees. In the steady state the P n (m) are related by the equations 
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P nP n (0) = f n (l . . . l)P n {n) = p(l - p), 
pin — m)P n (m) = pin — m + l)P n (m — 1) for m^0,n, 



which lead to the result 



P n (m) = (1 — p)/(n — m) for m < n, 



n—1 



P n {n) = l-(l-p) £ l/(7i -m) (5) 



m=0 

= l-(l-/>)£l/m- 

m=l 

The size distribution of forest clusters is 

n( S ) = P s+2 (01...10) = ^±^ 

= 7 T-(I-P)^ 2 . (6) 

(s + l)(s + 2) V HJ y 1 

This is a power law with the critical exponent r = 2. Eq. (§) leads to Y^f=\ n i s ) — (1 — p)/2, 
in agreement with Eq. (||). The size distribution of fires is oc sn(s) oc s -1 . The size 
distribution n e (s) of clusters of empty sites is 

n e ( S ) = P s+2 (10 ... 01) = = a(a + 1)(fl + 2) , (7) 



which is also derived in [12 



There is a characteristic length s max where the power law n(s) oc s~ 2 breaks down. 
During one timestep, a cluster of size s grows with probability 2p to a bigger cluster, and 
it is destroyed by lightning with probability fs. When the cluster size approaches s oc p/ f, 
the cluster is struck by lightning with a finite probability fs. This simple argument leads 
to s max oc p/ f and A = 1, as did the scaling theory in [0. The remaining critical exponents 
are p. = v = 1, as already derived in J7|. The following microscopic derivation shows that 
the relation for s max acquires a logarithmic correction factor. We calculate s max from the 
condition that a string of size n < s max is not struck by lightning until all trees are grown. 
When a string of size n is completely empty at time t — 0, it will be occupied by n trees 

6 



after T(n) = (1/p) J2m=i l/ m — ln(n)/p timesteps on an average. The mean number of trees 
after t timesteps is m(t) = n[l — exp(— pt)]. The probability that lightning strikes a string 
of size n before all trees are grown is fJ2t=^ m (t) — (f /p)n(hi(n) — 1) ~ (//p)nln(n). We 
conclude 

s max ln(s max ) ocp/f for large pj f. (8) 

Next we calculate the relation between the mean forest density p and the parameter / /p. 
The probability P n (n) in Eq. ([5]) cannot be less than zero. Since Ylm=i^/ m — hi(n) for 
large values of n, the mean forest density must approach the value one at the critical point 
via 1 — p oc 1/ ln(s max ). A more precise result for p is obtained from 

oo 

p = E sn ( s ) 

3=1 

Smax g CO 

= (1 " " } S (. + !)(• + 2) 

~ (1 - p) ln(s max ) + / sn(s)ds. 

-5 max 

With the scaling ansatz n(s) = (1 — p)s~ 2 C(s/s max ), the second term in the last line reduces 
to (1 — p) multiplied by a constant factor, and thus 

^ P = ln(s max ) + const. ~ \a(p/f) + const. 

~ \n(p/ f) for large p/f. (9) 

The forest density approaches the value one at the critical point. This is not surprising since 
there exists no infinitely large cluster in a one dimensional system as long as the forest is not 
completely dense. Combining Eqs. @ and @, we obtain the final result for the cluster-size 
distribution near the critical point 

n(s) ~ — — for s < s max (10) 

1 ; ( S + l)(s + 2)ln Smax max 1 ; 

with s max given by Eq. (g). 

Our computer simulations confirm the analytic results. The simulations were performed 

using the same method as in @. In Fig. | the size distribution sn(s) of clusters struck 



by lightning is plotted. The smooth line is the analytic result derived from Eq. (^j). It 
fits the simulation perfectly in the region s < s max ~ 1000. The bump in the cluster size 
distribution has the following explanation: consider the dynamics of a forest cluster. As 
long as the cluster size is s < s max , its average size increases each timestep by 2p(s + 2) (can 
be derived from Eq. (^)), i.e. its growth speed is s oc s. A tree spends a portion oc 1/s of 
its lifetime in clusters of size s, and consequently the growth speed is related to the cluster 
size distribution by sn(s) oc 1/s, which leads again to the power law n(s) oc s~ 2 . When the 
cluster size becomes comparable to s max , the growth speed increases no longer proportionally 
to s but slower than s since the largest neighbors of the cluster have already been destroyed 
by lightning. Therefore s 2 n(s) increases with increasing s for s > s max , giving rise to the 
bump in the cluster size distribution. When the forest cluster is so large that it is struck by 
lightning with a nonvanishing probability, n(s) decreases exponentially. 

The ratio ln(p//)(l — p)/p is given in Table | for different values of p/f. It approaches 
the value 1 with increasing p/f, as predicted by theory (see Eq. (|9|), the deviations from 
the value 1 for smaller p/f are due to the constant contribution in Eq. (0)). The ratio 
(p//)/(s max ln(s max )) is also shown in Table | for different values of p/f. Within numerical 
accuracy these ratios are identical, in agreement with our theoretical result Eq. (g). 

Finally, we calculate the temporal correlation function of the number of burning sites. 
The mean number N s (t) of trees that burn t timesteps after a cluster of size s is struck by 
lightning is N s (t) ~ 2(1 —t/s)9(s — t). The correlation function is then 

poo poo 

G(t)(x dsn(s)s dt N a (t)N a (t + t) 
Jl Jo 

oc (2s max / ln(s max ) - 3r(l - ln(r)/ ln(s max ))). (11) 

The Fourier transform of G(r) is G{uj) oc oj~ 2 {1 + const. In (cjs max )) for small to > l/s max . 
There is a nontrivial deviation from the uj~ 2 dependence in direction of l//-noise. Our 
simulations confirm the power law dependence (Fig. 0) but cannot discriminate between a 
uj~ 2 law and the complete expression Eq. (|TT|). 
We conclude with three comments. 
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(i) In the limit f/p—*0 only an infinitesimal portion of all forest clusters, all empty sites, 
and even all trees are not described by our results Eq. (|5])-(|7|) which therefore become 
exact at the critical point. 

(ii) The self-organized critical forest-fire model is in one dimension much simpler than in 
higher dimensions. The critical exponents are classical in one dimension, and the tree 
distribution on a string of size n < s max is stochastic. In higher dimensions, there are 
always trees left when a fire passes through a region. Therefore the tree distribution 
is not stochastic and the exponents are nontrivial. 

(iii) The only conditions for obtaining the result Eq. (||) are that trees grow stochastically 
and that a fire does not pass through a string of size n < s max as long as not all its trees 
are grown. The P n and consequently the cluster size distributions remain unchanged 
when the rules for the occurrence of lightning are changed in one of the following ways: 
(a) Lightning strikes all forest clusters of size > s max . This model is critical in the 
limit s max — > oo. (b) Lightning strikes every T timesteps the largest forest cluster on 
each interval of I sites. This version is critical for I — > oo with fixed T. (c) Lightning 
strikes a forest cluster of size s with a probability ef(s) which does not decrease with 
increasing s. This model is critical for e — > 0. While the size distribution of the forest 
clusters and of the clusters of empty sites remains unchanged in all these versions of 
the forest-fire model, the size distribution of fires varies considerably and is in general 
no power law. 
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TABLES 

TABLE I. Numerical results for the ratios ln(p//)(l — p)/p and (p//)/(s max lns max ) for differ- 
ent values oip/f. 



p/f 


3125 


6250 


12500 


25000 


50000 


Hp/f)(i-p)/ P 


1.13 


1.10 


1.06 


1.04 


1.02 


(p//)/(SmaxlnS max ) 


3.65 


3.61 


3.56 


3.62 


3.55 
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FIGURES 

FIG. 1. Size distribution of the fires for f/p = 2/50000 and L = 2 20 . The smooth line is the 
theoretical result which is valid for cluster sizes < s max . 

FIG. 2. The power spectrum G(u) of the fire density. The simulation parameters are 
f/p = 8/50000 and L = 2 20 . The straight line (theory) has the slope -2. 
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